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It is argued that the Brans-Dicke theory may explain the present accelerated expansion of the 
universe without resorting to a cosmological constant or quintessence matter. 



I. INTRODUCTION 

o ■ 
o ■ 

For several decades the general belief in cosmology had been in favour of a presently matter-dominated universe, 
C\| ! expanding as a oc i 2 / 3 , where a is the scale factor of the Robertson- Walker metric and t is the cosmic time. The 
• corresponding decelerated expansion (qo — —ad/a 2 = 1/2) was more or less compatible with all the cosmological 
tests. The problems with the standard cosmology were confined to the early stages of the evolution of the universe 
^| i and they were expected to be taken care of by the inflationary paradigm. But the recent observations regarding 
the luminosity-redshift relation of the type la supernovae up to about z = 1 revealed that the universe is in fact 
expanding at a faster rate, even possibly with an acceleration (qq < 0) jjj. This observation naturally leads to the 
search for the matter fields, hitherto unknown or neglected, which could introduce such a non-decelerated expansion. 
This matter is called a quintessence matter ( Q-matter for short). The list of possible candidates, being explored as 
^ ■ Q-matter, consists of a cosmological constant or a time dependent cosmological term, a scalar field with a potential 
giving rise to a negative pressure (2j at the present epoch, a dissipative fluid yielding an effective negative stress || or 
more exotic matter like a frustrated network of non-Abelian cosmic strings or a frustrated network of domain walls 
pj. Amongst the scalar fields chosen as the Q-matter, the "tracker" field rolling down its potential || appears to be 
the most attractive. But most of these investigations are effective only for a spatially flat (k = 0) model and thus 
the fine tuning or the flatness problem remains an unsolved one. The exception is the recent work by Chimento et 
al. Q , where a combination of a quintessence scalar field with a potential and a fluid with a bulk viscous stress has 
been shown to work for k = — 1. This approach also solves the coincidence problem in the sense that the ratio of the 
^-i density parameters for the normal matter and the scalar field asymptotically becomes a constant. 

The aim of this paper is to show that Brans-Dicke's scalar tensor theory can potentially solve the quintessence 
problem as it can lead to non-decelerating solutions for the scale factor for the present universe dominated by cold 
matter with negligible pressure. We would like to stress that this kind of solutions can be obtained without invoking 
a Q-matter or even dissipative processes. For a spatially flat model (k — 0) one can indeed have an accelerated 
, expansion, and even for a non-flat (k ^ 0) model, one can at least have a non-decelerated expansion. Very recently 
Bertolami and Martins || presented an accelerating model for the spatially flat (fc = 0) in a modified Brans-Dicke 
theory using a potential which is a function of the Brans-Dicke scalar field itself. In this work we show that this 
solution also can be obtained in Brans-Dicke theory without the potential. However, there are problems in bridging 
this result with the radiation-dominated decelerated universe for the same values of the Brans-Dicke parameter lu. In 
the next section we show that the solution to this problem lies in the natural generalization of the theory by allowing 
to be a function of the Brans-Dicke scalar field. 



II. FIELD EQUATIONS AND SOLUTIONS 

For a Friedmann-Robertson-Walker spacetime, with scale factor a(t), spatial curvature index fc, and assuming that 
the only matter-field is a perfect fluid, the gravitational field equations in Brans-Dicke theory are 
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a? + k _ p f _ a<t> w 4> 2 
a2 6 a ^ + 2 02' W 

a d 2 + fc = -p/ u> ft a4> <j> 

a a 2 2 <f> 2 a4> 0' 1 j 

where w is the Brans-Dicke's constant parameter, while pf and pf are the density and hydrostatic pressure respectively 
of the fluid distribution. These two latter quantities are connected by a barotropic equation of state Pf = (if — 
7/ being the (constant) adiabatic index of the fluid. 

The wave equation for the Brans-Dicke scalar field reads 

0+^ = ^-^, (3) 
a Zlu + 6 

and this combined with the field equations ( [l]) and ( |^) leads to the matter conservation equation 

Pf + 3^(pf+p f ) = Q. (4) 

Assuming that at the present epoch the universe is filled with cold matter of negligible pressure (i.e., dust), we put 
Pf = in the last equation and obtain the relation 

Pf=Poa~ 3 , (5) 
where po is a constant of integration. With pf — and pf being given by equation ( ||), the wave equation becomes 

(« 3 #- P ° 



2uj + 3 7 

which immediately yields the first integral 

a 3^ = ^o^. (6) 

We are primarily interested in a power law non-decelerating solution for the scale factor in keeping with the recent 
observations. So we take 

a = a t a (a > 1) (7) 



where clq is a positive-definite constant. 

Equation (^) can now be integrated for <j> as 



*2-3q 



(8) 



v ag(2w + 3)(2-3a)' 

Inspection of the field equations ( Q) and ( ||) reveals that when k ^ 0, the only possible solution in this simple 
power function form is that with a, = 1, i.e., 

a = a t, (9) 

and 

Po 



ag(2u + 3)i' 



(10) 



So one has a scenario where the present deceleration parameter vanishes (q^ = 0), i.e., today the Universe is in a state 
of uniform expansion. Although this very simple model does not yield an accelerating universe for a non-zero k, it 
expands with no deceleration either and may sufficiently fulfil the requirements set by the recent observations on the 
distant supernovae [Q. 

From (Q) the consistency condition (for a = 1, i.e., qo = 0) 
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w = -2(l + 4) (11) 
a 

follows. The latter equation indicates that at least for k — or k — +1 , u is negative. This is perfectly consistent 
with equation ( |l0|) which requires that w < — 3/2 as is a positive quantity. 

If the universe is spatially flat (k = 0), this model yields several possibilities including accelerated expansions for 
the universe. With k = 0, the equation ( ^J), put back in the field equations yields for dust (pf = 0) the general 
solution as M a = 2(ui + l)/(3w + 4), i.e., 



and the scalar field takes the form 



= aoi 2(«+l)/(3*+4) j (12) 

(3cj + 4) p ^ 2/(3 " +4) 



2ag(2a; + 3) 



(13) 



Now it is easily seen that for different negative values of u>, one can generate different accelerating solutions. Such as, 
for u> — —5/3, we get a = 4/3 and (f> cx i -2 , which is the solution presented by Bertolami and Martins Q. In this 
case, go = —1/4. For ui — —8/5, a becomes 3/2 and <p cx £~ 5 / 2 , giving go = — 1/3, i.e., an accelaration rate higher 
than the Bertolami-Martins' solution. For different values of u> in the range —2 < w < —3/2, this model yields a 
host of accelerating solutions for a spatially flat universe. In fact this negative value of uj explains the behaviour of 
the model. A sufficiently negative u may effectively lead to a negative pressure (see equation (|^)) and thus drive a 
positive acceleration or uniform expansion without the violation of the energy condition by normal matter. 

The solutions and the corresponding value of w can be used in equation ( to calculate the age to of the uni- 
verse. It turns out that the age is of the order of Hq 1 . For the uniformly expanding spatially flat universe, to is 
exactly equal to H^ 1 ~ 1.62 x 10 10 yr. The higher the acceleration, the older the universe is. For go = — 1/4, 
*o = f #0 

As Brans-Dicke theory is a varying G theory, the rate of this variation has been checked in this model and it appears 
to be compatible with the observational limit. For instance, with go = — 1/4, (G/G)o = 3Hq/2 < 2.5 x 10~ 10 which 
is safely below the upper limit of 4 x 10~ 10 yr -1 set by observations (see || and references therein). 

It is important to note that although for the some negative values of ui, the model can lead to decelerating expansions 
for the Universe when it is radiation-dominated, in that case the value of ui would be less negative than —3/2, i.e., 
—3/2 < to < 0. The indicated range of values of ui, i.e., —2 < ui < —3/2, which drives an accelerated expansion 
for a matter-dominated model, does not produce a consistent model with the radiation-dominated epoch. Thus the 
suggested model seems to badly spoil the big-bang nucleosynthesis scenario. One way out of this problem is to start 
with a modified version of Brans-Dicke theory where the parameter ui is a function of the scalar field <f> rather than 
a constant In this case equations ([j]), (||) and (Q) remain in place, but the wave equation for the scalar field (||) 
has an additional term and becomes 

3ai = Pf -3p f 

Y a 2w + 3 2lu + 3 v ' 

With the equation of state for radiation (pf = \pf), equation ( |l4| ) immediately gives a first integral 

^(2c + 3) 1 / 2 = 4, (15) 

where A is an integration constant. It turns out that there are functional forms of ui for which the universe expands 
with a deceleration such as like a cx <2 in this case so that the primordial nucleosynthesis can be successfully explained. 
A simple choice like 

2cj + 3=(0-1) 2 (16) 
when used in equation ( |l5| ) along with a cx i 1 / 2 yields the integral 

(0-l) 2 = ^, (17) 
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B being a constant of integration. Equations and ( jl7| ) clearly shows that u decreases with time towards a 
constant value —3/2, which indeed produces an accelerated expansion for a late dust-dominated universe as we have 
seen. 

In general, a choice of u) as a polynomial in <£> is 

n 

2w + 3 = 5^i4 i ^, (18) 

i=l 

where Uj > and Aj's are constants appears to solve this problem. A varying lu theory with a nonminimal coupling 
between gravity and the scalar field different from the Brans-Dicke theory (4> 2 R as opposed to <fiR in BDT) may also 
explain the late time behaviour of the universe as pointed out very recently by Bartolo and Pietroni jllj . It is also 
worthwhile to note that Sen and Seshadri [ ft2| have shown that with small negative values of u>, it is possible to obtain 
growing modes for the density perturbations for a universe with a power law expansion rate. 



III. A CONFORMALLY TRANSFORMED VERSION 



In Brans-Dicke theory, the folklore is that uj is positive, or at least larger than —3/2. There are two reasons for this 
belief. The first reason is that in the weak field limit, the Newtonian constant of gravitation G is given by, 



G = 



2uj + 3\ 1 
2lo + 



and for —2 < u < —3/2, G becomes negative. But it should be emphasized that this is only in a weak field limit. So 
in the model presented in the present work, G does not become negative, only this weak field approximation does not 
hold. In the full non linear Brans-Dicke theory, the effective Gravitational constant is actually given by 

which is indeed positive as (f> is positive. It deserves mention that for the accelerating solution obtained by Bertolami 
and Martins || also, u> is actually negative. 

The second reason for the belief that ui > —3/2 is that the energy related to this scalar field is proportional to 
(2ui + 3). To see this, it is better to effect a conformal transformation 

9n» = (f>9nv 

In the transformed version, G becomes a constant, but one has to sacrifice the equivalence principle as the rest mass of 
a test particle becomes a function of the scalar field |l3| . So the geodesic equations are no longer valid and naturally 
the physical significance of different quantites are indeed questionable in this version of the theory. Nevertheless the 
equations in this version gives an insight in comparing the energies of different components of matter. The equations 
( [!])-( ||) in the new frame look like 

3(a 2 + k) _ 2w + 3. 2 

=2 = Pf + — a — ^ , (19) 



„ ^2 

2a a + k 2ui + 3 



a 

and the wave equation ( [|) transforms to 



= -pf — (20) 



^ + 3^ = ^/, (21) 
a 2lo + 3 

where an overhead bar indicates quantities in the new frame, ip — lncf>, and a 2 — 4>a 2 . The density and the pressure 
of the normal matter in this version are related to those in the original version in a simple way as pf = <ft~ 2 Pf and 
pf = <j)~ 2 pf. It is clearly seen that the contribution to the energy density by the scalar field is given by 
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2^ + 3 ; 2 

Pt = b 2 - (22) 

Here t/j looks like a massless minimally coupled scalar field and at least formally behaves as a perfect fluid with a stiff 
equation state, p^ = p^. (Note that the corresponding adiabatic index is 7^ = 2). Equation ( |22| ) explains why u> is 
usually taken to be larger than the —3/2 from the consideration of the positivity of energy. 

Our intention is to get a sufficiently negative pressure from some source so that we can get go < 0, and here we 
achieve that by means of a scalar field </> which has a negative energy. But in our case ^ is a geometric field unlike 
the non-gravitational fields such as normal matter or the quintessence field, and so its having a negative energy is 
not pathological in that sense. We rather save the normal matter from the unappealing feature of a negative pressure 
and yet get a non positive-definite deceleration parameter. 

This conformally transformed version of the theory gives an insight regarding the solution of the flatness problem 
too. The density parameter is defined as f2 = p/3H 2 , where 

2lu + 3 • , 

P = Pf + P4> = Pf + — 4 — V> > 

is the total energy density. The subscripts / and 4> refer to the normal fluid and the scalar field, respectively. H = a/a 
is the Hubble parameter in the new frame. 
The Bianchi identity yields the relation 

p + 3jHp = 0, (23) 
where 7 is the average barotropic index defined by the relation 

7^ = 7/0/ +7^0^, (24) 

where 0, = (if + fi^. 



Using equations ( |19|) and ( 23) one can write 

n = n(n-i)(37-2) J ff. (25) 

This equation shows that (l — 1 is indeed a solution, but I d(l/d(l j must be negative at Cl = 1 if this (l value is to 

be a stable solution of ( |25|). This requires that 7 should be less than 2/3. As the ratio of the density to pressure of 
the fluid remains the same in the two versions of the theory, the adiabatic index 7/ also remains the same. Because 
of our choice of pf — and p<f, — p ( f >1 we have 7/ = 1 and 70 = 2. So from equation ( |24|) one has 

7= %±|i. ( 26) 
As 2uo + 3 is negative, Q$ is also negative and one can achieve 7 < 2/3 provided 

(if < 4 1 o 1 . 

So fl = 1 is indeed a stable solution in this model depending on the relative magnitudes of the energies of matter and 
the Brans-Dicke's scalar field. With (l = 1, we have = —k/a 2 = and thus the flatness problem can be solved. 

It is true that although the conformally transformed version (popularly known as the Einstein frame) is appealing 
from the point of view of the computational simplicity, the quantities lose their physical significance for reasons stated 
at the beginning of this section. But it must be emphasized that the character of k remains the same and if it is zero 
in this version, it is so in the physical original version of the theory as well. For a spatially flat universe (k — 0) we 
can thus construct accelerating models for the Universe and reproduce the Bertolami-Martins solution as a special 
case and we do not have to invoke an additional self-interaction term for this. Even for the non-flat cases, we can at 
least obtain a non-decelerating model with qo = 0, whereas almost all the quintessence models produce solutions only 
for a spatially flat universe with the only exception of the dissipative fluid model recently discussed by Chimento et 
al. ||, which works for open universes as well. For a discussion of the flatness problem in Brans-Dicke theory see the 
paper by Levine and Freese p3. 
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IV. CONCLUDING REMARKS 



Brans-Dicke theory proved to be extremely useful in solving some of the outstanding problems in the inflationary 
universe scenario with the possibility of an "extended inflation" |l5f| . And now although the theory might relinquish 
some of its natural attraction as it fails to produce general relativity as an infinite lo limit |l6| ] contrary to the earlier 
belief, once again it appears to be able to account for some outstanding features of our present Universe like the 
accelerated expansion and the coincidence problem by its own right without having to invoke dissipative processes or 
exotic fields. 

It should not be overlooked, however, that lo has to pick up a low negative value (i.e., — u> ~ 0(1)) in order to 
solve the acceleration and coincidence problems by a single stroke (e.g., qo = —1/4 implies a oc t 4 ' 3 and consequently 
lo = —5/3). This squarely conflicts with the lower limit imposed on |cj| by solar system experiments, namely |u>| > 500 
jl7j . Nevertheless in "extended inflation" the model of La and Steinhardt worked provided that lo takes a value close 
to 20 [ fl5| which is also embarrasingly lower than the mentioned astronomical limit. Therefore it remains a problem 
to find a suitable compromise between astronomical observations and cosmological requirements. Another problem 
is that the constant negative value of lo does not produce a consistent radiation model which explains the primordial 
nucleosynthesis. We have shown that a varying lo theory like Nordtverdt's jl0| can give rise to a decelerating radiation 
model where the big-bang nucleosynthesis scenario is not adversedly affected and lo evolves to the small and negative 
constant values required for the late time acceleration in the matter-dominated epoch. A thorough survey of varying 
lo theories is perhaps warranted which may bridge the decelerating radiation universe with the accelerating matter 
universe while the local inhomogeneities might locally give rise to high values of lo to be consistent with astronomical 
experiments. 

However, attempts to solve the aforesaid problems outside Brans-Dicke theory do not fare much better. The 
introduction of a cosmological constant of the order of the critical density lacks of a solid support from quantum field 
theory flgfl , and arguments based on the anthropic principle are anything but convincing fl9| . On their part models 
based on quintessence suffer from the problem of unwanted long range forces and that they cannot be as homogeneous 
as it should [§0) . 
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